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ON OBLIQUE GROWTH OF TREES UNDER THE 
ACTION OF WINDS 


I. OPATOWSKI 
THE UNIVERSITY OF CHICAGO 


It is shown that an inclined regular growth of tree trunks under 
the action of prevailing winds may be explained as a response of the 
plant to the mechanical action of the wind, in conformity with the 
Schwendener theory of shape of plants based on the concept of maximum 
strength. A possible mechanism of such response is suggested. The con- 
dition of maximum strength determines the inclination of the tree trunk 
on the basis of purely mechanical considerations. A numerical example 
computed for a palm agrees with observational data. 


Since the fundamental investigations of S. Schwendener (1874) 
(cf. also Holtermann, 1909) a very large amount of experimental ma- 
terial has been published on the relations between mechanical proper- 
ties of organs of plants and their form and functions; the response 
of plants to different types of mechanical stimuli has been also inves- 
tigated (for a review of these topics see W. Schwarz, 1930, 1931 and 
W. Rasdorsky, 1931, 1937). The ordinary geotropism is a well known 
type of such response, and other types have been observed directly or 
induced artificially through various forms of mechanical stress. A 
tree trunk with an elliptic cross-section having its major axis in the 
direction of prevailing winds or an increased growth of the trunk and 
of the roots in the regions which have to resist higher loads (Bern- 
beck, 1913; Kroll, 1913; Knight, 1803; Preda, 1901; Rasdorsky, 1925) ; 
changes in anatomic structure induced by bending, the so-called kamp- 
totropism (Biicher, 1906; Neubert, 1911), exhibiting, for instance, 
a thickening of cell walls or an eccentric formation of annual rings 
(Hartig, 1901, pp. 51-57), are all examples of responses of plants to 
stimuli consisting of mechanical stresses (cf. Bower, et. al. 1939, 
chapter x; Haberlandt, 1924, chapter iv; Ungerer, 1926, pp. 247-251). 
The biophysical mechanism of such responses may be connected per- 
haps with the pressure gradient, as it arises, for instance, in bending of 
a branch through an increase of compressive stresses in one part of the 
branch and occurrence of tensile stresses in the remaining part (cf. 
Bernbeck, 1924) ; the importance of pressure differences for the growth 
of a tree is suggested by the fact that about 30 to 60 per cent of 
green wood consists of liquid substances and consequently in the 
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presence of a pressure gradient, a substantial intercellular flow may 
occur with a possible change of relative concentrations within sin- 
gle cells, if the permeability of the cell walls exhibits a specific dif- 
ferentiation for the various constituents of the cell. 

According to S. Schwendener, the shape of a plant organ having 
a mechanical function is such as to achieve a maximum of strength. 
This idea, which was known already to Galileo, has been discussed 
many times since Schwendener. It is the purpose of the present note 
to examine this principle of maximum strength in connection with 
the action of winds on trees. 

The mechanical effect of winds on the trunks of trees has sub- 
stantially two different forms. One occurs when an unusually strong 
wind bends the tree into a plastic range, so that a very large portion 
of the tree trunk may acquire permanently an abnormal position, the 
trunk continues, however, to grow vertically. A beautiful example 
of such effect is given by L. Preda (1901) in an old juniper having 
about one-half of its trunk completely horizontal. The same type of 
effect is shown by the famous “praying” palm of Faridpore in India. 
This is a date palm whose trunk was bent through 60° by a storm and 
whose crown exhibits diurnal movements from an erected position 
in the morning to a lowered position in the evening, the latter effect 
being a response to environmental changes of physical character 
(Bose, 1929, chapter xxi). Many other interesting examples of ab- 
normal shapes of trees due to a plastic bending of the trunk may be 
found in the publications of F. C. Schiibeler on the flora of Norway. 

A different type of wind effect consists of an oblique growth 
of the whole tree trunk under the action of prevailing winds. A 
beautiful photograph of a group of palm trees growing in this way 
was published recently in the National Geographic Magazine (Moore, 
1945). [For other examples see e. g. F. C. Schiibeler (1885)]. The 
question arises in connection with this second type of effect whether 
such regular oblique shape of a trunk should be explained also on 
the basis of a plastic deformation which has occurred at a certain 
moment during the life of the tree; or whether under a continuous 
action of the winds a biological response of the tree could affect its 
growth in such a manner as to give an oblique shape to the tree. If 
the latter is the case, we would conclude from Schwendener’s theory, 
by limiting ourselves to the mechanical action of the wind alone, that 
the inclination of the tree should correspond to a minimum stress. 
We examine the physico-mechanical characteristics of the tree and 
its environment in which such explanation is possible. 

With reference to the already quoted group of palms we consider 
a cylinder of length H as an approximate geometric model of the 
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trunk. Let the horizontal cross-sections of the cylinder be circles of 
radius R, and the angle of inclination of the cylinder with respect 
to the vertical axis be v. The wind velocity is variable with the 
height, particularly near the ground; mathematical formulae express- 
ing this variation are known (s. e. g. Sverdrup, 1939; Brunt, 1941, 
pp. 251-252) and could be used here; for our approximate computa- 
tions, however, we assume an average velocity V constant along the 
entire height H. Theoretical considerations and experimental re- 
sults indicate that the tangential component with which a wind may 
act on a circular cylinder is negligible and that the normal component 
WwW: is 
W:=Ncos*(h+ v), 


where h is the angle between the direction of the wind and the hori- 
zontal, h being taken positive if the wind is directed upwards and 
negative if it is directed downwards; N is the force of the wind when 
it blows perpendicularly to the cylinder (cf. Finsterwalder, 1902, pp. 
164-166; Feld, 1945) and is given by the formula (cf. Betz, 1925, pp. 
370-373) : 

N=RHV, 


where the coefficient k, varies from 0.004 to 0.006 for H/R varying 
from 20 to ~,N being in kg ,V in km/h, R and H in meters. We will 
use, however, a much higher value of k; to take into account the 
roughness of the trunk. Let C be the weight of the crown acting at 
a horizontal distance c from the top of the trunk. We assume that 
the force with which the wind acts on the crown is representable by 
a formula of a well-known type in aerodynamics; viz.: 


W.—k,.AV? cos?h , 


where A is the area of the vertical projection of the crown, k, is a 
numerical coefficient depending on the shape of the crown and the 
force W, is assumed to act horizontally at a vertical distance K from 
the top of the trunk. Substantially the same type of formula has been 
used by P. Jaccard (1919) in his studies on trees. Calling the weight 
of the trunk T , the compressive stress arising at the base of the trunk 
is due to the compressive force 


r=—TsW,sin? +C, 
and to the bending moment around the center of the base: 
M= (H/2)(W,;+Tsinv) + (¢+ Hsinv)C + (K+ Heosv) We. 
The maximum compressive stress S at the base of the trunk is, ac- 
cording to known formulae of the strength of materials, 


S= (F + 4MR")/(a B*), 


44 OBLIQUE GROWTH OF TREES 


so that by the previous formulae 
aR?S/H =k,RV2(2H + Rsin v) cos?(h + v) + 2(T + 2C) sin v 
(1) 
+4W,cosv+ H7[(T+C)R+4(cC + KW.)]. 


Assuming that the wind force acting on the crown is independent of 
the inclination v of the trunk, a minimum value of S, if it exists, 
occurs when v satisfies the equation 


d(x R? S H-* V2) /dv =k,R[ F cos v cos? (h + v) 
— (2H + Rsin v)sin2(h+ v)] (2) 
+ 2V(T + 2C) cosv—4k.A cos?*hsinv=0. 


By its physical meaning, the Schwendener principle of maximum 
strength requires that, as v increases from zero, the function S(v) 
has a minimum as its first extreme. Therefore, S(v) must be a de- 
creasing function at v = 0, and from equation (2) we obtain the 
condition 


(cosh — 8l sinh) cosh + 2[(T + 2C)/(k:R?V?)] =0, 


where | = H/(2R) is the slenderness of the trunk. We obtain a ne- 
cessary condition for the fulfillment of the above inequality by chang- 
ing in it 2C into C, but 


So= (T+ C)/(a R?). 


is the stress S of the same tree when growing vertically in absence 
of winds. Consequently, we obtain the condition 


(cosh — 8 l sinh) cosh + 22 So/(k:V?) <0. (3) 


From this inequality we conclude that a sufficient slenderness of the 
tree, low compressive stresses in the trunk, and high velocities of up- 
wards directed winds, are necessary conditions for the ability of the 
plant to control its growth under the action of winds in conformity 
with Schwendener’s principle. In the case of downwards directed 
winds, the above inequality cannot be satisfied because all its terms 
are positive. Experiments carried out in Germany by O. Lilienthal 
in various conditions of terrain up to a height of 10 meters show 
that the average winds are directed upwards with an inclination of 
about 34° to the horizontal (Betz, 1925, p. 367). With such small 
value of h, normal trees cannot satisfy the inequality in question. 
In exceptional terrains, however, as, for instance, in mountains or 
on sea beaches, the direction of the wind may be sufficiently steep or 
its velocity sufficiently high to make a biological response of the tree 
in the sense of Schwendener’s principle possible. 


I, OPATOWSKI 45 


The condition (3) may be written also in the form 


sin(2h—a) > [1+ (42S,)/(k,V?)] (1 + 64/)+, (4) 
where 
sina = (1 + 647) = 1/(8I), (5) 


l being very large with respect to 1. If we consider the maximum 
compressive strength of the wood o, and introduce the safety factor 
of the tree fy = o./S,, corresponding to vertical growth in absence 
of winds, we have (cf. Opatowski, 1944, pp. 113-114) 


So = 473 X 10‘G/f, kg/m? , (6) 


where G is the specific gravity of oven dry wood based on the volume 
in the green state. The right hand side of inequality (4) must be 
< 1; this gives a condition independent of the direction of the wind; 


lf,k.V? > 74 X 10°G, (7) 


which is obtained with the help of equations (5) and (6). The bio- 
physical meaning of conditions (3) and (7) is the same. The aero- 
dynamic coefficient k; depends on the shape of the tree and on the 
roughness of the trunk; it is probably not smaller than 0.02. The 
slenderness / of the trunk may exceed 100. The safety factor f, must 
be high, because it is referred to the action of the weight of the tree 
alone, and the tree has to resist also strong winds: f, may well exceed 
100 or 200. The specific gravity G varies between about 0.3 and 0.7. 
On the basis of these numerical data we may see from relation (7) 
that an oblique growth of a tree in the sense of Schwendener’s prin- 
ciple of maximum strength is possible only in exceptional conditions 
involving winds of about 50 or 75 km/h at least. The requirement 
of such exceptionality of conditions is in agreement with common 
observation. It should be pointed out, however, that our calculations 
are based on the assumption of a cylindrical tree trunk; a trunk in 
the shape of a truncated cone, or having some shape of maximum 
strength could lead to slightly different numerical conclusions. It 
should be noted also that the conditions of slenderness and of high 
safety factor are contrasting with each other. For this reason, it is 
worthwhile to eliminate from relation (7) the product lf, by means 
of the radius of the trunk R. This can be done in the following man- 
ner: putting p = C/T and calling m the ratio of the weight of the 
. green wood to its weight when oven dry, we easily obtain from equa- 
tion (6) . 
m(1+ p)R1If, = 2365. (8) 


Consequently, from relation (7) 
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m(1 + p)GR < 82 X 10°k,V?, (9) 


where, as before, R is in meters, and V in km/h, and the humidity 
factor m varies between about 1.3 and 2.5. The meaning of inequality 
(9) is that the effect under consideration may be observed only on 
thin trees which do not have heavy crowns. 

We proceed now to a calculation of the inclination of a tree, on 
the basis of equations (1) and (2). Consider, for instance, a palm 
about the size as represented in the already quoted article of W. R. 
Moore: R = 0.075 m and H = 15 m, which corresponds to a slen- 
derness of 1 = 100, p = 0.1, A = 3 m?; assume for the specific grav- 
ity G and for the humidity coefficient m the values given by L. J. 
Markwardat and T. R. C. Wilson (1935) for a Florida palm: G = 0.37, 
m = 2.34. This corresponds to a weight of the trunk and of the crown 
of T = 229.5 kg, C = 22.95 kg. Equation (8) gives the safety factor 
corresponding to a vertical growth of the tree in the absence of winds: 
fp = 122.5. Putting k, = 0.04, the minimum value of V for which the 
condition (9) is satisfied is slightly less than 75 km/h. We assume 
V = 80 km/h, although much higher wind velocities have been also ob- 
served (cf. Feld, 1945). Equations (4) and (5) give the minimum 
value of wind inclination for which the effect may be observed: h is 
somewhat above 30°. We take A = 38°. Assuming an aerodynamic 
coefficient for the crown of k. = 0.0005, which is about ten times 
smaller than for a circular plate, we obtain the force of the wind act- 
ing on the crown W, = 5.96 kg. With these values the equation (2) 
for the inclination of the tree derived from Schwendener’s principle 
of maximum strength is 


cos v cos? (38° + v) — (400 + sin v) sin (76° + 2v) 
+ 382.5 cos v — 16.56 sin v=0, 


whose numerical solution gives v = 214°. This value is probably 
slightly smaller than the inclination of the palms in the photographs 
of W. R. Moore. However, we have neglected the vielding of the 
ground. 

It is convenient now to calculate the actual safety factor of the 
tree, in these conditions, and to see how much it gains through the 
inclined growth. The compressive strength of wood is much higher 
in pure bending than in pure compression. Formulae of L. J. Mark- 
wardt and T. R. C. Wilson give the compressive strength in compres- 
sion o, and in bending o, (cf. Opatowski, 1944): 


a, — 473 X 10°G,, op = 1237 X 10*G*/4 © (10) 


in kg/m?. Equation (1) gives the actual compressive stress S: the 
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terms of the right hand side of equation (1) due to pure compression 
are 
k.R?V? sin v cos?(h + v) + H7?R(T+C), 


the remaining terms being due to pure bending. We put the compres- 
sive stress at the base of the trunk in the form 


S=S,.+S,, (11) 


where S, and S, represent respectively the parts due to pure com- 
pression and to pure bending, both under the combined action of the 
weight of the tree and of the wind. We take c = K = 0, that is, we 
assume that the weight of the crown and the resultant wind force 
acting on the crown pass through the top of the trunk (cf. the photo- 
graphs of W. R. Moore) ; these forces do not affect sensibly the re- 
sult in any way. We obtain from equation (1), using our numerical 
values, 


=iSe-— 14° S, = 404.8 forv=0, 
; (12) 
S-=— 1.4; S; = 269.6 for v = 214° 


the units being kg/em*. It is seen that S, is negligible with respect 
to S;, so that in calculating the strength the second formula (10) 
can be used and the safety factor is f = «/S. Using equations (11) 
and (12), we obtain f = 0.88 for v = 0, that is, assuming our tree 
is vertical, and f = 1.32 for an inclined growth of the same tree ac- 
cording to Schwendener’s principle of maximum strength. The 50% 
of reduction of stress achieved through an appropriately inclined 
growth is here of a substantial importance for the tree, since it can- 
not live with f < 1. The safety factor of 1.32 is not very high, it is 
calculated, however, for a strong wind, and it should also be kept in 
mind that our theory is based on a cylindrical model of a tree, which 
is very inefficient from the viewpoint of strength. A model with a 
suitably variable cross-section would yield a higher value of f. Such 
a model could be determined by combining the ideas of this paper 
with a theory developed by M. Esser (1946) which gives a shape of 
a trunk of maximum strength. The latter, however, should be extend- 
ed to include the action of the wind on the trunk. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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THE UNIVERSITY OF CHICAGO 


Following a previous paper, equations are derived for the most 
probable time of firing of an efferent neuron in terms of the intensity E 
of excitation of the afferent pathway, when E£ is either constant or any 
given function of time. The equations are not differential equations, but 
in integral form. It is suggested that ©, correspondingly j, represent 
the number of excitatory, correspondingly inhibitory, terminal bulbs ex- 
cited within the period of latent addition at a given most probable time. 
The relation between the suggested theory and the old one, based on dif- 
ferential equations for © and 7 is discussed. 


In a previous paper (Rashevsky, 1945b) we suggested an inter- 
pretation of the « and 7 factor as the number of excitatory or in- 
hibitory terminal bulbs, which is most likely to be excited simultane- 
ously at a given time t. We suggested a non-linear differential equa- 
tion which reduces as a very rough first approximation to the stand- 
ard equation for e and 7. The validity of this equation is, however, 
rather dubious from the point of view of theory of probability. 

In the present paper we shall derive some more rigorous prob- 
abilistic equations. We shall find that a closed finite expression for 
_the excitatory, correspondingly inhibitory, factor is thus obtained 
directly, without any recourse to differential equations. For those 
_ used to handling the problems of the mathematical theory of the cen- 
tral nervous system in the standard way, based on the differential 
equations for e and 7, this may appear strange. But we must remem- 
ber that after all a differential equation is only a mathematical ex- 
pedient, which possesses some heuristic value. Neither in physics 
nor in any other mathematical natural science do we ever test ex- 
perimentally a differential equation directly, but only some finite ex- 
pressions derived from the differential equation. In many cases, as 
in diffusion problems, eletrodynamics, hydrodynamics, the differen- 
tial equation has the advantage of a greater generality. But in some 
cases, as, for instance, in the case of central excitation, we really do 
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not gain anything mathematically by using the differential equation 


Of eae tyee (1) 
dt 


‘= 
instead of its general integral 


t 
e(t) see &—A { E(t) e dt h (2) 
0 WW iktaee ea 
Both give us « as a function of ¢ for any prescribed E(t). 

Consider the same set-up as before (Rashevsky, 1945b), namely, 
an efferent neuron whose threshold in terms of the minimum number 
of terminal bulbs which must be excited is n,. Let N afferent fibers 
synapse with this neuron, each fiber having only one terminal bulb 
on the neuron. Let, < N. Let the average frequency » of discharge 
be the same in all afferents and let 6 be the period of latent addition. 

The probability P(n) for any n fibers to bring their impulses to 
the efferent neuron within the interval 6 is then given (Rashevsky, 
1945b) by 

N! 


(N—n) In!" 


We shall call the firing of any terminal bulbs during the in- 
terval 6 “event n’’. 

The number N varies with the intensity E of the afferent exci- 
tation, and therefore for a varying E the number N is a function 
of time ¢. 

The passing of each interval of time 6 may be considered as a 
trial, and the probability of the event » to occur during a given trial 
is then given by equation (8). During the 7-th consecutive trial, that 
is, in the 7-th interval 6 after the beginning of the afferent excitation, 
the number WN has the value N; , and hence P(”) is a function of both 
i and ”, which may be written P;(n). 

The probability of the event » to occur once and only once dur- 
ing m trials, when its probability to occur during the i-th trial is 
equal to p; , is given by 


p(1,m) =p, (1— pz) «++ (1— Dm) 


P(n) = (6 ¥)" (3) 


4 
sf (1 — p:) 2 (1 — ps) -+-(1 — Dm) Hees, (4) 
The r-th product of the above expression may be written 
D+ i=m : 
II (1— pi). (5) 


L— De ie 
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Consider now the logarithm of the product in expression (5). 
We have 


log I (1— pi) =S log(1—pi). (6) 
i=A1 4=1 
The interval 6 is very small (about 0.5 ms), and m is therefore 
usually very large. If t denotes the time corresponding to a given 
fixed value of m, then 
t 
t=m 0; 6=—. (7) 
m 
In our case 
pi—P;(n). (8) 
If, as is usually the case, P;(n) = p, << 1, we may expand the 
logarithms in expression (6), preserving only linear terms, and thus 
obtain 


log I (1—p:) =¥ log(1—p:) =— Spi. (9) 


For very large values of m, the above may be approximated by 
t 
= f p(r)d7r, (10) 
where t denotes the fixed time considered, and 7, the integration vari- 
able, is the variable time, of which p(r) is now a function. We have 
p(r)dr=pi. (11) 
Hence 


li O—p) =exp -f roa: | (12) 


i=1 
Now the fraction p,/(1—,) in expression (5) becomes simply 
p(r)dz+. Hence every term of the sum in expression (4) is now equal 
to 


| exe|— J paz ]| pla)d; Seek 
and therefore, writing p(t)for p(1,m): 
p(t) =| exp [-f/poe-]] free. (14) 


For p(r) =const.=p, this reduces to 
p(t) =pte*’. (15) 
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The quantity p(t) has a maximum for 


==. (16) 
p 
or, using equations (8) and (11) and putting p= p; ,dr+=0: 
= (17) 
P;(n) 


which is identical with equation (2) of the previous paper (Rashev- 
sky, 1945b). 

If the condition p; < < 1 is not satisfied, we must evaluate the 
sum (4) directly. For the case p; = Const. = p,, equation (4) re- 
duces to: 


p(1,m) =mp.(1— 7p)”. (18) 
This has a maximum for 
1 
m = — ——————_. (19) 
log (1 — Do) 
For very small values of » , equation (19) reduces to 
1 
Po 


which, because of relations (7) and (11), is identical with expres- 
sion (16). 

Now let us ask for what value of ¢ the probability p(t) has a 
maximum in the general case, given by equation (14). We have 


ee = ptt exp [- f pera] 
—| J pear| p(t) exp ie | bene) (21) 

The requirement 
ar (22) 


gives 


r(t) | ie J vca:} exp |- fome| =0, (23) 
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or. 


if p(r)dr=1. (24) 


For p(r) = Const. = p, this reduces to equation (16). 

The number N is a function of the intensity H(7) of the afferent 
excitation. Hence N is a given function N(r) of the time. Because 
of equations (8), (8), and (11), together with dr = 6, the quantity 
p(r) is thus a known function of the variables + and n; thus p(r) = 
(7, %). Hence for every value of n , equation (24) gives us the most 
probable time for which n bulbs will be once and only once excited all 
within the interval 6. Putting n =, , we find the most probable time 
for the efferent neuron to fire once and only once. 

For a constant FE, the value of N is constant, and the most prob- 
able time of firing of the efferent neuron is given by 


} ri) (N —7,) !n,! 


P(n,) (6»)™. WN! 


(25) 


In actual experimental determination of reaction times, the ex- 
periment is repeated many times under constant conditions, but each 
observation yields a different value of the reaction time. All those 
values are scattered around a most probable one, which is considered 
as the reaction time. We may therefore consider equation (25) 
as giving the reaction time as a function of the intensity of the stimu- 
lus S in a simple reflex are. For large values of N, and forN >> n,, 
that is, for stimulations highly exceeding the threshold, ¢, varies ap- 
proximately as N-™. Hence, ifn «E « S, then t, « S™. Actually, 
it is found (Landahl, 1939, Rashevsky, 1940, chap. ix) that 4, « S7 
and sometimes even t, « (log S)-1. Since most likely n, >> 1, the ex- 
pression (25) leads to a too rapid decrease of t, with increasing S , 
than is compatible with experiment, if N « S. But actually, as S in- 
creases N will tend to a limiting value N,, the total available number 
of fibers in the pathway (Rashevsky, 1938, chap. xxv). The N—S curve 
will rise for large values of S slower and slower. The actual law of 
variation of N with S depends on the distribution function of the 
thresholds. Empirically we may try to put 


as 
b+S 


It would be of interest to try how well reaction time data, ana- 
lysed by H. D. Landahl (1939; Rashevsky, 1940, chap ix) from the 
point of view of the differential equations « and 7, will agree with 


(26) 
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equations (25) and (26) for a proper choice of the constants a, b, 
and n,. 

The present theory gives us not only the reaction time as a func- 
tion of the stimulus intensity, but also the frequency of the observed 
deviations from the most probable value. For a constant stimulus 
equation (15) gives us directly the probability of a response for any 
time ¢. 

Equation (24) defines an implicit function n(t). Putting 


Fin,t) = [ reoa= fi o@.mar, (27) 
we have 
oF 
dn _ at 
A (28) 
an 


From equation (27) we have, putting 


Ty ee Sat 
aK apa ee (29) 
oF oF . 
—S— t : ——— a T,n dr. : 
aot ms = foe (7,2) (30) 
Hence 
dn o(t, n) 
4 hoe Ls SF (31) 
f dn (7, n) dz 
or 
ae o(t, 2) 
n (7, 2) dr = — ———_., 
i, gn (7,2) ee (32) 
dt 


Taking partial derivatives of both sides with respect to ¢ and rear- 
ranging, we find 


dn : dn \? : dn 
6(t,n) 4, (t,m)( 5) (tn) =0. (38) 


This is the differential equation determining the number of ter- 
minal bulbs firing within an interval 6 as a function of the most prob- 
able time of that event. For practical purposes the integral expres- 
sion (24) is much simpler. 
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The number n may be identified with «, correspondingly 7. If 
equations (25) and (26) can represent intensity-time relations in 
reaction times for constant stimuli, then the present theory is empiri- 
cally equivalent to the standard one, because no actual quantitative 
comparison of theory and experiment for variable stimuli has been 
made. To what extent the two theories may be equivalent mathe- 
matically, can be decided only by a further mathematical study. If we 
set the right side of equation (2) equal to h, for a given E(t) it de- 
termines the firing time ¢, of the efferent neuron. We may compare 
the roots of this equation with that of equation (24) forn =n, . 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara Abbott Memorial Fund of the University of Chicago. 
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Let S denote the intensity of a stimulus, v the frequency of dis- 
charge of a fiber, and », the integral frequency of discharge of a nerve 


pathway. The relation between », and S can be derived when the relation 


between » and S, as well as the distribution function of the fiber thres- 
holds are given. It is investigated under what conditions the functions 
v7(S) and »(s) will be both of the same exponential form, and it is shown 


that this can happen only if S exceeds the highest threshold of the path- 
way. 


Peripheral nerve fibers are excited by applying a stimulus of 
sufficient strength to their receptor organs. For rapidly adapting 
fibers, a few discharges will occur, and then no further excitation 
will occur unless the stimulus is increased. For slowly adapting fibers, 
a steady volley of discharges at a frequency y» sec will be evoked for 
some length of time. Touch receptors are an example of rapid adap- 
tation, while muscle receptors such as the flower spray receptor organ 
are instances of slow adaptation (Fulton, 1943). We propose here 
to consider pathways of parallel fibers which are so slowly adapting 
that they may be approximately considered as non-adapting, all sub- 
ject to the same stimulus S. 

For these fibers, the frequency discharge y will be given by some 
functional relation 


y=9(S,h), (1) 
where h, the threshold of the fiber, is the smallest stimulus capable 
of producing excitation, and the function ¢ is monotonic increasing 
with S but may have an upper limit. In addition we have 

6(S,h) =0, Sih, (2) 
If we have N, fibers in all in the pathway, we define the total 
frequency of the pathway 7 as 


=D %- 3 yea (3) 
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N. Rashevsky (1945) suggests that vp is the measure of the total ex- 
citation E of a pathway. If h’ is the lowest threshold in the pathway 
and N(h), h 2 h' is the distribution function of thresholds in the 
pathway, then assuming that ¢(S,h) is the same function for all 
the fibers in the pathway, we have 


Ss 
(8) = f o(S,h)N(h)dh. (4) 
BA 
Making the assumption that ¢(S,h) is given by 
il 
a(S, h) =7 lesen], CoSeen 


o(S,h)=0, S<h 


a function suggested by N. Rashevsky (1940) as possibly applying 
to some nerve fibers, equation (4) becomes 


at : e-2es : 
r(S) == | N(h)dh—— ee" N(h) dh. (6) 
h* h’ 


It is of interest to see what conditions would be imposed of neces- 
sity on the distribution function N(h) in order that »7(S) have for 
all values of S the same form as the function of equation (5) as- 
sumed to hold for an individual fiber, i.e., let »r(S) be given by 


yp (S) =k[1 — eS"). (7) 
Denoting by M(x) and L(a) any two functions satisfying 
dM (x) dL (x) 
—_=N : — ever : 
ie (x) , de N(x) ; (8) 


and combining equations (6), (7), and (8), we obtain 
1 
e[L— e6*] == (M(S) —M(h) —e***[L(S) ~L(h’)]}. (9) 
Differentiating both sides of equation (9) with respect to S gives 


bke-b(s-h') — ae-°8 | 7, (S) ir L(h')] ’ (10) 
hence 


bk 
L(S) — L(h’) =— e’ eae-d)s | (11) 
a 
Differentiating both sides of equation (11) with respect to S gives 


bk 
eves N(S) meget (a6 — b) e2°S ets (12) 
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from which we obtain, as a necessary condition on the distribution 


function N(h) in order that »; be a function of the same form as the 
function ¢ of equation (5), the following: 


bk 
N (h) =— (a6—b) ec?" , (13) 
a 


Now if the function N(h) of equation (13) is introduced into equa- 
tion (6), the result is a function »7(S), given by 


k 
yr (S) eae [ (ag — b) + be) — age] , (14) 

a 
This function has a different form from the function of equation (7), 
hence we must conclude that there exists no distribution function 
N(h) which for all values of S will make »7(S) a function of the form 


1 
k{1 cannes ere ii when ca) — ol. =e foe Mah 2 


Any nerve pathway has a maximum threshold h,, , i.e, N(h) =0 
forh >h,. Then for S > h,,, equation (6) takes the form 
v(S) =A — Bes , (15) 


where A and B are given by the equations 


1 = 
= if N(h)dh =N./6 


ya (16) 
B=- i) em" N(h) dh 
Ba bc 
If we denote by p the quantity 
1 B 
=— lo oad (17) 
£ aé a 
then equation (15) may be written as 
N 
A penicee erckaaasgis (18) 


a function of the same form as the function of equation (5). 

Thus for a stimulus S great enough to excite all the nerve fibers 
of the pathway, the total frequency of discharge vr of the pathway 
will depend on S in the same manner as the discharge frequency of 
an individual fiber. This result is independent of the distribution 


function N(h). 
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If hm is very small, equation (18) may be simplified for values 
of S reasonably close to hy by expanding e-*°*”) and dropping all but 
the linear terms. We thus obtain 


yp(S) =N.a(S—p). (19) 
Applying the mean value theorem for integrals to the second equa- 
tion (16), we obtain 


B= even (20) 


where h is between h’ and h,,. From equations (20) and (16) we 
have 


B = gash (21) 
A > 
hence 
nlp 2 (22) 
—— lor —=— pp. 
ag ~ A a 


Thus p is between h’ and hy, . 
If we assume that »;(S) is a linear function S, 1.e., 


-ye(S) = KS+F, (23) 


introducing equation (23) into equation (6) and making use of equa- 
tion (8) we have 


KS+ P=" [M(S) —M(h’)] —= 


eS 


[L(S)—L(h’)]. (24) 
Differentiating both sides of equation (24) gives, after rearranging 
terms, 
’ K ‘ 
L(S) —L(h) ane (25) 
Differentiating both sides of equation (25) gives 


N(S)=6K, (26) 


as a necessary condition on the distribution function N(h) in order 
that vp (S) be linear. On the other hand, if N(h) of equation (26) is 
introduced into equation (5), we obtain for »7(S) 


yp(S) = K(S—h’) = [1 —- see 7) . (27) 


This function is certainly not linear and at most is approximately 
zero for values of S close to h’. Thus we may conclude that »7(S ) will 


; 
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not be linear for a very large range of S for any function N(h). 

Much work in the mathematical biophysics of the central ner- 
ous system is based on the assumption that for the ranges of stimuli 
considered, »7 is approximately a linear function of S. We see from 
preceding results that for this assumption to be valid, it is necessary 
to assume that the maximum threshold h,, of the nerve pathway is 
very small and that the quantity a@(S—p) be sufficiently small. If 
yr(p) is close to zero, then equation (18) and the additional condi- 
tion that vy = 0 for S < p corresponds very well to the variation of 
yp with S. If the distribution function N(h) were skewed to the 
right considerably, then the value of »7 would be small and would not 
rise rapidly until h approached very close to p. In such a case as 
this, we could use equation (19) with good results. Since 6 is a meas- 
ure of the refracting time for nerve fibers in the pathway, we would 
obtain better results with nerves of very low refractory time. 

This paper was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of the University of Chicago. 
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The theoretical shape of the trunk of a tree is determined by pos- 
tulating that the maximum compressive stress due to the action of grav- 
ity and wind on an inclined tree is constant for all cross-sections of the 
trunk. In the general case, the shape is obtained by successive approxi- 
mations. In some particular cases, the shape is given explicitly. 


Introduction. A theoretical study of the shape of tree trunks 
can be based on various postulates as to its determining cause. It 
can be based, for instance, on 


I. Considerations of sap conduction, where the trunk is con- 
sidered mainly as a conveyor of metabolic substances. 

II. Mechanical considerations, where the trunk is considered 
mainly as a support of the crown. 

III. Considerations of growth, in which the rate of increase in 
diameter of the trunk is considered more important than the actual 
diameter at any time. 


We cite a few names: N. Rashevsky (1948a, b) studied I and II. 
S. Schwendener (1909) and I. Opatowski (1944a, b; 1945a) studied 
II, and P. Jaccard (1919) studied a combination of I and III. 

Besides the three factors mentioned above, many other causes 
influence the shape of trunks. An ideal treatment of the problem of 
shape would take all factors into consideration, but for the simplicity 
of ‘the mathematics involved, only one or two factors should be con- 
sidered at one time. The results so obtained should then be compared 
with actual data. Experimental data are not yet precise enough to in- 
dicate the relative importance of the leading different theories. 

In the present and in a subsequent paper, the viewpoint II will 
be used. 

For our study, trees consist of a trunk and a crown, the trunk 
being the part of the tree located between the ground and the base 
of its lowest branch, and the crown being the remainder of the aerial 
part of the tree. When a primary branch is considered, this branch 
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will also be separated into two parts, namely, the portion contained 
between the trunk and the first secondary branch, and the remainder. 
These two parts will be called the stem and the crown of the branch 
respectively. We will use the language corresponding to trunks, but 
it must be kept in mind that our considerations apply also to primary 
branches because these branches can be considered as very inclined 
trees attached to the trunk instead of being rooted to the ground. 

Trunks are considered to be homogeneous solids with circular 
symmetry. The whole tree is supposed to have a plane of symmetry. 
The action of gravity and of wind is considered. The force of wind 
on the crown is supposed horizontal and in the plane of symmetry 
of the tree. 

By “cross-section 2” we mean a cross-section of the trunk by a 
plane normal to the axis of the trunk, and at a distance z from the 
base of the trunk. We determine the shape of the trunk by postulat- 
ing that the maximum compressive stress o(z) on each cross-section 
equals some constant o: 

o(z) =c. (1) . 


Sheering stresses are generally small in comparison to compres- 
sive stresses and therefore will not be considered. The assumption 
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(1) could be justified, for instance, as follows: 

The shape of a theoretical trunk is such that it can support a pre- 
assigned crown with the greatest economy of material. Mechanical 
function or biological activities impose an upper bound o on the stress 
allowable in the trunk. The greatest economy of material is obtained 
by giving a minimum radius to each cross-section. To a decrease of 
radius corresponds an increase of stress, and the minimum radius is 
obtained by setting the maximum stress o(z) on each cross-section 
equal to the upperbound oc. 

Instead of assuming a greatest economy of material, equation 
(1) could also have been justified by assuming a maximum height of 
the tree, or a greatest safety factor against breakage, or by showing 
how stress affects growth. 


Calculation of the stress. We use the following notations (cf. 
Figure 1): : y) 


= 
=i 
aa! 


0 = point of intersection of axis of trunk with ground 
l= length of trunk 
G= force of gravity on crown 
W = force of wind on crown 
Y = horizontal distance from 0 to the line of application of G 
Z = vertical distance from 0 to the line of application of W 
z= distance from a variable point of the stem to the cross- 
section through 0 
r(z) =radius of cross-section 
6= angle of inclination of trunk with the vertical 
w= specific gravity of stem 
o(z) = compressive stress at the lowest point of the cross-sec- 
tion at z. 


In order to obtain always the maximum stress on the lowest 
point of each cross-section, we suppose that 


620 
Gceosé—Wsin62=0 (2) 
G(Y —Isin 6) + W(Z—lIcos 6) 20. 


We consider a force F in the vertical plane containing the axis 
of the trunk and applied to a point of the crown or of the portion of 
trunk contained between the crown and the cross-section correspond- 
ing to a certain z. Denote by F, the projection of F’ upon the axis 
and by M(F) the moment of F' with respect to the center of the cross- 
section. Then the compressive stress due to F' at the lowest point of 
the cross-section z is given by (Timoshenko, 1940, p. 226) 
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FB; 4M (F) 
a(z) = = 


a r2(z) are(z) 


When several forces act on the tree, the resulting stress is the 
sum of the stresses due to the different forces. 

We get for the stresses oc, or, ow, due to the gravity of the 
crown, to the gravity of the trunk, and to the wind respectively, the 
following values: 


Geosé 4(Y—zsin6)G 


bison nr? (Zz) nr (Zz) 
cos d.@. fi r(C)de Asin6.w. ['(C—z)r(C)de 
on (2) = + —__—_————_, 
os (2) (2) 
—Wsinéd 4(Z—zcosé)W 
ow (z) =— SS o 


769? (2) n 7? (z) 


The total stress is 
o(z) = o¢(Z) aia or (2) =e ow(2) 
G @ fi 72 (£) at 
= cos 6 2 aerate at an 
m1 (z) 7? (z) 


4 2G 4 ow f! a | 


uA (3) 
L. 2% 7*(@) r*(z) 


E G [ sin 6 ee | 
Y | - +W | ——— + -—————_ | 
zr? (z) nevis) ne 7a) 

The last member consists of four terms which consider respec- 
tively the gravity for a perfectly vertical tree, the influence of inclina- 
tion of the trunk, the influence of eccentricity of the crown and the 
action of the wind. 

A numerical example. The purpose of this example is to show 
the relative importance of the different contributions to the stress 
o(z), and also to see to what extent our postulate that «(z) is inde- 
pendent of z, see equation (1), is justified. The units used in the 
calculations. are meters and tons. 

A complete set of data on any one tree has not been found. 

Some data of geometrical nature on the “epicea no. 1 de Hasliwald” 
are given by P. Jaccard (1919, pp. 3-5). The missing data will be 
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obtained by more or less speculative methods. The height of the tree 
is 38 meters and the height of the trunk is 1] = 15 meters. The weight 
of the trunk is calculated by assuming a specific gravity w = 1 and 
the weight of the crown is obtained by assuming o(3) = o(15) in 
the absence of wind and inclination. The force of the wind, W = 3.5 
tons, is obtained from an estimation made by P. Jaccard for an in- 
tensity of wind such that it exerts a force of 50 kg on a plate of one 
square meter. This is a strong wind of a little less than 100 kilometers 
per hour. The point of application of W has been taken at the lower 
third of the crown (Z = 22), although Jaccard for some unknown 
reason, takes it still lower. We get then 


TABLE I 


a “ men 
E é E Ey s 3y8 
8 3 Sn Pais) =: ZOE 
re = oe 23 os 9 HF 
wm i oo P+ Bo y Bou 
a ® Se mo HE ww Bo HS 
om ood QE Ooo & as Bas 
sea ESE Y oSS 6 OR) 35 
i og Bas Ey a x oa 
8 || 3 & || ws 36” o% 248 
wn 
3 z 2 Bs Es gue 
rn a = 23 ne tal 
s S) a 
| | | 
15 466 | 170 | 0 | 18.5 | 1250 
| | | | | 
13 | 495 192 | 854 | 182 = 1320 
| 
11 | 510 | 204 | -748 | 19.0 1480 
| | | 
9 | 524 | 215 | 1.160 | 20.0 | 1610 
| | 
if | 554 | .240 | 1.610 | 19.8 | 1580 
| | | {| 
5 | .588 | .270 | 2.100 | 19.4 | 1500 
| | a 
3 | 635 | 315 | 2.680 | 18.5 | 1330 
| | | i] 
1 | 75 | 440) | 


3.390 14.8 | 900 


Calculated weight of the crown: 


weight of crown _ cross-section at 15 G 170 


ot) 


weight of tree ~ eross-section at3’ G+2.68 .315 
G= 3.14 Tons. 
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We will now compare the order of magnitude of the four terms 
of the last member of equation (3). The first term is of the order of 
20 tons per square meter, as shown in the preceding table. We will 
calculate for what inclination, eccentricity, and wind intensity the 
other terms of the last member of equation (3) are also of the order 
of 20. 

To calculate 6, we have the equation 


4 2G 4a yh (ft —2z)7r?(C) dt 
sine ~- ff |=». 


a 1 (2) (2) 
A rough estimate of the brackets will be obtained by taking in the 
first term z equal to zero, taking as an average value 7(z) = .30 me- 


ters, and taking the integral equal to 


i I 
l 
«[ froe= <a [ roe 
1 1 


l 
ss X mass of the trunk = 5 X 3.89 = 17 meter tons. 


Co|—™ 


We have then 


Brawl? , 
sin 6 =20; sind=.0075, @—degree. 
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The eccentricity Y is given by the equation 


i 4 G 4 3.14 
20=Y- = Y -— = # —— . 
EET are 150 Y, as 138 meters 


The numbers in the last column of the preceding table are about 80 
times larger than those of the preceding column. To get numbers of 
the same order of magnitude in both columns, we would have to take 
a wind of intensity 50/80 = .& kilograms per square meter. This 
wind has a velocity of about 12 kilometers per hour. 

Conclusions from the numerical example. We see in the table 
that the stress due to gravity is fairly constant and therefore justifies 
the postulate o; = constant [ef. equation (1)]. The stress due to 
wind has a greater variation. The postulate o = constant, is there- 
fore less justified. The quantities og and oy have the same defect of 

being too large in the middle of the trunk. Therefore, a postulate of 
the form og + k ow = constant, where & is a positive constant, and 
which takes into account a combined action of gravity and wind, is 
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not more justified by our example than the postulate og = constant. 
The example shows also that for inclinations greater than half 
a degree, eccentricities greater than 13 centimeters and wind veloc- 
ities greater than 12 kilometers per hour, the stresses due to inclina- 
tion, to eccentricity, and to wind are each greater than the stress in 
a perfectly vertical tree in a calm atmosphere. It is thus evident that 
none of all these factors can be neglected. 

To show the importance of eccentricity and inclination, we cite 

a formula found on pages 230-231 of S. Timoshenko (1940): 


d 
o=o.| 1+ 4° | ; 
fp 


which gives the maximum stress o on an arbitrary circular cross- 
section of the trunk in terms of the stress o, in a perfectly symmetric 
tree of same weight and same radius of cross-section, of the radius 
ry of the cross-section, and the horizontal distance d of the center of 
gravity of the tree to the center of the cross-section. 

Another mechanical factor which may be important but will not 
be considered, is the elastic stability of the tree. 

Theoretical form. A combination of expressions (1) and (3) 
gives 


o 7? (Z) =@ cos 67 (z) [ e@a Se 


Ll 
4 sin f (¢ —z)r2(t)dt + (4) 
G cos 6— W sin 6 G sin 6+ W cos é GY+ WZ 
eS (ge A 
Jt 7U LA 


This integral equation can be solved by the successive approximations 


19 (z) = arbitrary function ; 


I 
o Tani (Z) = @ COS 6 1n(Z) I 72,(C) dae + 


dosing [ (—2)r%(C)at (5) 


G cos 6— W sin 6 G sin é@+ Wcosé@ GY + WZ 
 —— — ——————_ z + 4 --——_; 
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The validity of this method of solution can be proved in the usual 
way, as it is done, for instance, by E. Goursat (1929, chap. xix). 

If for all z we have 7,(z) = 17:(z), then the sequence 7,(z) is 
non-decreasing. To show this, we remark that from 7,(z) = 7,(z) and 
from inequality (2) it follows that the second member of the mid- 
dle equation of system (5) is smaller for n = 0 than forn=1. A 
comparison of the first members gives then 7,(z) = 72(z). Proceed- 
ing inductively, we find that 7,(z) = tiu(z) for each n, which is 
what we wanted to show. The inequality 7,.(z) = 7,(z) is satisfied 
in particular by 7.(z) = 0 and by the functions 7,(z) which we will 
consider later. 

Physical considerations suggest a set of approximations a little 
different from those given by expressions (5). The approximation 
Tnui(z) is calculated so as to obtain a constant stress when we take 
as weight of the trunk not its actual weight, but the value determined 
by the approximation 1r,,(z). The equations of these approximations 
are obtained from expressions (5) by replacing in the second mem- 
ber of the middle equation the first r,(z) and the last r,(z) by 1nu (2), 
while the 7,(¢) under the two integral signs are unchanged. This ap- 
proximation method involves the solution of an algebraic equation of 
third degree, but does not involve the solution of an integral equation 
or differential equation. 


Solutions for particular cases. The equation (4) will be solved 
exactly in some particular cases. 

I. w=0O. In other words, the weight of the trunk is negligible. 
The integral equation (4) reduces then to an algebraic equation of 
third degree. If in addition we set 0 = G = 0, we get 


4W 
r (z) mey (2s8)s (6) 


Thus when the action of wind alone is considered, our trunk is shaped 
as a cubic paraboloid with concavity downward. 

Il. 6=0,Y=0,W=0. In other words, we have a perfectly 
vertical tree in the absence of wind. Equation (4) reduces to 


eso f m(eyae + <, (7) 


By differentiation we get 
dr? (z) 
gels Lod 
dz 


=— or (z) ’ 
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and we obtain as a solution of equation (7) 


v2 (z) Sheen (8) 
to 
This exponential form has been considered by many authors and is 
sometimes in fairly good agreement with observations. 
We obtain also a simple second approximation if we take 6 = 0, 
Y = 0, W = 0 for our first approximation and only 6 = 0 for our 
second approximation. We get then, when 7,(z) is given by equation 
(8), 


4 4G 
eae) Wa (Ze) ey 
wo ao 


e 3/2 : ; AW 4AGY 
el Bele ese(!-2)/20 4+ —__ (FY — z) + 


Ito ao eo 


We see that this approximation 7*,(z) is the sum of the exact values 
of r*(z) obtained by considering the action of gravity alone, as in 
equation (8), of wind alone, as in equation (6), and of eccentricity 
alone. 


Hit. .¢ = me W = 0. In other words, we are studying a hori- 


zontal branch in the absence of wind. Then equation (4) reduces to 


. 4G 
om(eyase f (6 — 2) 0) ae ee CY — 2) (9) 


An exact solution of equation (9), involving hyperelliptic inte- 
grals, will be discussed elsewhere. The solution of equation (9) is 
simple in the following two cases: 

a) When w = 0, that is, when the weight of the crown of the 
branch is large nalattvely to the weight of the Sune, primary 
branch. Then 


4G 
r3(z) —=— (Y-—2). (10) 
ao 
b) When, on the contrary, G = 0, that is, when our primary 


branch is denuded of secondary fwoniclicn: then equation (9) has the 
solution 


Fe) =—* ea (11) 


as can be verified by substitution. 
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It should be noted that such a branch not only has no biological 
utility, but, as remarked by I. Opatowski (1945b, sec. 5), physically 
cannot exist. A calculation of its deflection leads to infinite values. 
This objection can be made whenever at the tip r(z) is an infinitesi- 
mal of at least second order in (/—z). 

If we take equation (11) as a first approximation and substitute 
into the second member of equation (9), we get as a second approxi- 
mation 


8 4 
(2) =| ~ (He ee (12) 


15¢ 


We see that this approximation 7*,(z) is the sum of the two exact 
solutions (10) and (11) obtained by considering respectively the 
weight of the crown alone and the weight of the branch alone. 

The author is indebted to Professor N. Rashevsky and Dr. I. 
Opatowski for their contributions in the preparation of this article. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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The activity of a neural net is represented in terms of a matrix 
vector equation with a normalizing operator in which the matrix repre- 
sents only the complete structure of the net, and the normalized vector- 
matrix product represents the activity of all the non-afferent neurons. 
The activity vectors are functions of a quantized time variable whose 
elements are zero (no activity) or one (activity). Certain properties of 
the structure matrix are discussed and the computational procedure 
which hele, from the matrix vector equation is illustrated by a specific 
example. 


Let a neuron be assumed to have the following properties. Its 
activity is of the “all or none” character. Unless it is an afferent 
neuron it can be brought into activity only by the action of other neu- 
rons which have terminal bulbs upon it, provided these act upon it all 
within the period of latent addition. Its activity lasts throughout its 
‘time of action,” beginning at the end of a synaptic delay. The syn- 
aptic delay is measured from the time at which sufficient conditions 
have first been met. More specifically, we shall assume that the syn- 
aptic delay time is the same for all neurons, and that compared with 
this, the period of latent addition, the time of action and the conduc- 
tion times are negligible. For convenience this interval will be taken 
as the unit of time, ard we shall treat time as being quantized. 

By a synapse shall be meant the relationship between the termi- 
nal bulb of one neuron and a dendrite or cell body of another neuron. 
If such a synapse can activate the second neuron under any condi- 
tion, it shall be referred to as an excitatory synapse and the terminal 
bulb as an excitatory bulb. If a synapse can ever prevent activity of 
the neuron acted upon, it shall be referred to as an inhibitory synapse 
with an inhibitory “bulb.’”’ We shall assume that one inhibitory syn- 
apse can act to neutralize the effect of one excitatory bulb in the same 
synapse group. Evidently total inhibition can be obtained by n* in- 
hibitory synapses where n*is a number equal to or greater than the 
total number of excitatory synapses in the group of synapses acting 
upon a common neuron. 

Let 6, be the least number of excitatory terminal bulbs acting 
simultaneously upon neuron N;, which is necessary for the neuron 
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to act a unit of time later. 

Let E*; be the number of excitatory bulbs from neuron N; ter- 
minating upon neuron N;, , and J*; the number of inhibitory bulbs, and 
let BY; = E*, — I*;. 

Define 
fix = BY;/O « (1) 


Evidently f;, may be positive, negative, or zero. It is zero if 
E*, = I*; or if N; does not synapse with k , in which case E*; = 1*; =0. 

Let the group of synapses implied in f;, for all values of 7 form 
a synapse group on neuron k. Assume that for every synapse group 
the effects of all terminal bulbs is simply additive. If a subgroup did 
not add with the rest of the group but only added within itself, then 
the neuron k may be replaced by three neurons. Two of these, N,- 
and N;-, are associated with the respective subgroups only and each 
acts on the third, which has a threshold less than the number of bulbs 
terminating on it either from N, or N;--. The third neuron acts 
upon the same set of neurons as does N;. The time delay introduced 
can generally be compensated for in the rest of the net. If there are 
several subgroups which overlap, this introduces a more complex 
reduction. We shall hereafter assume that each neuron has but one 
afferent synapse group. 

Let a net of n neurons be divided into afferents, efferents, and 
internal neurons. We shall only mean by an afferent a neuron not 
acted upon by any neuron in the net under consideration. Similarly 
an efferent is a neuron which does not act on any other neuron in the 
net. We may, however, refer to the afferents as receptors and the 
efferents as effectors. Let the p receptor neurons be N,, N2, --- N,, 
--- Np, the internals be Np. , --- Ni, --- Now and the « effector neu- 
ron be Nye, ++: Ne, +:: Nn. Define a structure matrix F as a square 
matrix having the row and column indices 1, 2, --- » corresponding 
to the m neurons such that 


Latin Ptl...4... (M-€41)...€...n 


(2) 
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Each element in row 7 determines which neurons are acted upon 
by neuron N; and in what manner. Similarly each element in column 
k determines which neurons act upon N; and in what manner. The 
matrices Fz, Fx, F;, and F;, appearing in equation (2), determine 
respectively the relationships receptor-internal, receptor-effector, in- 
ternal-internal, and internal-effector. Define the matrix R as the 
nm X n matrix, obtained from F by substituting zeros for all elements 
except those of F;. Define in a similar way the matrices X , I, and 
E. We shall assume that no neuron acts upon itself so that fx, = 0 
for all &, and F and F, have all diagonal elements equal to zero. 

If any row a in F contains only positive or zero elements, then 
Nz is a purely excitatory neuron. If any row # in F contains only 
negative or zero elements, then Nz is a purely inhibitory neuron. If 
both positive and negative elements occur in a given row, the corre- 
sponding neuron may be referred to as a mixed neuron. 


Define the (1 X ) row matrix or vector a, by 
a(t) = (@,,+++@,,+++Ap » Ups °°? Dj 9°99 Ae y An-e+u1 3*°*Ge,*** An), (3) 


where any element a; is 1 or 0 depending on whether N; does or does 
not act at the time t. The vector a(t) may be referred to as the ac- 
tivity vector at the time t. This vector may be written as the sum 
of three (1 X n) vector components, r,i, e, the receptor, internal, 
and effector components having the respective set of elements a,, 
a; , and a, only, and zeros elsewhere. The scalar quantity v;(¢) given 
by the sum 

V(t) => BY /h:., (4) 

jep 

taken over a class § of neurons which is defined as the class of all 
neurons synapsing on N; which are active at t—1, gives a measure 
of the net excitation affecting the neuron N;. The quantity v;(¢) 
may be expressed as 


mS fee Ss a1); (5) 


jep j=1 


because of the definitions of f;, and f. A vector v(t), whose compo- 
nents for k > p are the values of v;,(t) and for k = p—the values of 
a,,(t), can be expressed as 


v(t) =r(t) + a(t—1)F. (6) 


In order to normalize v(t), let G be a post-operator on a row vector, 
such that, if [v G], is the kth component of the vector [v G] (¢), 
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(vUGhk=1 if 21, 
7 
[vG],=—0 if. Meds 7) 


Since G operates only on a vector appearing just to the left of it, 
and not on a matrix, the expression (vF')G may be written uFG. The 
parenthesis could not be eliminated if a pre-operator were used. 

Equation (6) may now be written 


a(t) =v(t)G=[r(t)+ a(t— 1) F]G. (8) 


The post-operator G will not in general be distributive with respect 
to vector addition. It is not associative with v and F. However, in 
equation (8) the vectors a(t —1)F and r(t) satisfy [a(t —1)F]: X 
[r(t)], = 0 for all values of k, in which case it can be readily seen 
that G is distributive. Also rg = r. Thus we have the following 
theorem: 

The activity of any net represented by a structure matrix F ts 
determined from the afferent stimulation and its activity at the be- 
ginning of the prior interval of time according to the equation 


a(t)=r(t) +a(é—1)FG. (9) 


Equation (9) leads to the recursion formula 


a(t) =r(t) + {r(é—1) + [r(t—2) +--- [r(2) + [r) + 
(10) 
a(0)FG] FG] ---FG | FG} FG. 


The vector quantity between the first and last brackets in expres- 
sion (10) is simply the vector a(t — 1). 

Thus if the structure of a net is known together with a sequence 
of r’s, r(0), r(1),---r(+ — 1), and the initial activity of the inter- 
nal neurons i(0), it is possible from equation (10) to determine the 
activities of the net for any time t from 0 to 7, a(1), a(2), --- a(r). 

Because of the character of the structure matrix F, equation (9) 
may be written as a pair of equations 


i(t) =[r(t—1)R+i(t—1)G, (11) 
e(t) = [r(t—1)X + i(t—1)E]G, (12) 
from which one may determine successively i(t) and e(t) for every 


t. For both i(t) and e(t) formulas similar to expression (10) can 
be written. 


From equation (11) it is evident that the sequence of i’s can be 
determined from a knowledge of the afferent-internal structure and 
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internal-internal structures, together with the sequence of afferent 
activities and an initial internal activity pattern. On the other hand, 
to determine the sequence of e’s, that is, the pattern of the efferent 
activity, one must also know the rest of the structure of the net. 

Every row and column of F, excluding the first p columns and 
last « rows, contains at least one non-zero entry; otherwise, it repre- 
sents an afferent or efferent. If only one non-zero entry occurs in 
any column, it may be replaced by unity, for if it is less than one, the 
neuron of the corresponding column can never act, and thus this neu- 
ron should be deleted. Furthermore, there is no restriction to set an 
element equal to one, if this element is greater than one. 

For each term in the expansion of the determinant of F; for 
which a class of q non-zero factors have subscript pairs forming a 
permutation of order gq, there is a closed chain of q neurons. If each 
factor is = 1, it is a functional closed chain. Let D; represent any 
term in the expansion of the determinant of F’; in which zero elements 
of the term are omitted and the non-zero elements are written as 
fix. Then for every D; in which any ciass of the subscript pairs 
forms a permutation, there corresponds to this a closed chain of 
neurons, the number cf neurons being the same as the order of the 
subscripts. If M, is the p X » matrix obtained from the matrix IJ by 
deleting all rows and columns except those associated with these p 
neurons, and if the expansion of the determinant contains only one 
non-zero term, then the closed chain is simple if the subscript pairs 
form a permutation of order p , and no proper part of this permuta- 
tion is itself a permutation. Otherwise, each neuron is a member of 
a simple closed chain, each chain being disconnected with any other 
insofar as these p neurons are concerned. If each element f;; within 
a permutation is 2 1, the closed chain is functional. If any row or 
column in J contains more than one non-zero entry, and each contains 
at least one non-zero entry, there is at least.one complex closed chain 
or one pair of connected closed chains. If the (« Xv) matrix F; is a 
permutation matrix and if r(0) = r(1) --- = 0, then after a finite 
time the sequence of e’s is cyclic and of order « or some submultiple. 

To give an example of the application of the present method, we 
exhibit the matrix for the circuit illustrated in Figure le of the paper 
by W. S. McCulloch and W. Pitts (1948): 
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es 
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If the net is initially at rest, then according to W. S. McCulloch 
and W. Pitts the condition for the activity of neuron 3 at t = 0, is 
either that neuron 2 acts at t = —3 but not at t = —2; or that neuron 
1 acts at t = —1. That is, the sequence of r’s, writing only the 4, 
components, r(—3) = (0, 1), r(—2) = (0, 0), r(—1) = (0, 0) as 
well as the sequence r(—1) = (1, 0) is adequate to produce activity 
in neuron 3. We shall assume that r(t) =0 for allt 20. Ifr(—3) = 
(0, 1) then from equation (9) a(—2) = (0, 0;1,—1;0,})G= 
(0, 0; 1,70; 0; 0), a(—1) = (0, 0; 0, 1; 0, 0), and a(0) = 
(0, 0;0,0;1, 0), so that e(0), writing only the a. components is 
e(0) = (1, 0), that is, neuron 3 acts but 4 does not. Similarly if 
r(—1) = (1,0) then from equation (9) e(0) = (1,0). Note that for 
both sequences a(1) = 0. 

The condition for neuron 4 to act as t = 0 is that neuron 2 acts 
at ¢ = —2 and att = —1. If r(—2) = (0,1) and r(—1) = (0,1), 
then from equation (9) a(—1l) = (0, 1; 1,0; 0,0) and a(0) = 
(0,0; 1,0; 0, 1), so that neuron 4 acts but 3 does not. Note that 
for this sequence a(1) = (0,0;0,1;0,0) and a(2) = (0,0;0,0; 
1,0) so that neuron 3 always acts as a unit of time after discontinua- 
tion of continuous stimulation of neuron 2. 

In the above example unnecessary restrictions were imposed for 
simplicity. Evidently, if the net is initially at rest and x denotes a 
quantity which is either zero or one, then the sequence of r’s given by 
r(—3) = (#, 1), r(—2) = (2, 0), r(—1) = (2, x) or the sequence 
r(—1) = (1, «) results in the efferent activity vectors e(—2) = 
(x, 0), e(—1) = (0, 0) and e(0) = (1, 0) for the first sequence 
and e(0) = (1, #/2)G = (1, 0) for the second. On the other hand, 
the stimulus pattern r(—2) = (x, 1) and r(—1) = (0, 1) results in 
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re(—1) = (x, 0) and e(0) = (0,1). 
_ This work has been aided in part by a grant from the Dr. Wallace 
CC. and Clara A. oe Memorial Fund of the University of See 
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